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A microscopic picture of the "preparation" of a crystal 
to the transition to liquid state at the approach to melting 
temperature is proposed. Basing on simple crystallogeometric 
considerations and the analysis of the computational results 
for corresponding anharmonic characteristics an evaluation is 
given for the magnitude of atomic displacements leading to 
the appearance of close packed Bernal pseudonuclei in the 
solid phase. A physical meaning of the classical Lindeman 
criterion and the mechanisms of the formation of free volume 
in the crystal phase are discussed. 
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Despite a great amount of cdWsfied experimental and 
theoretical information (see, e.g.,ErEl) a classical question 
what is the melting is still unsolved. The melting temper- 
ature Tm itself is determined from thermodynamic con- 
siderations and, e.g., for alkali metals it can be calcu- 
lated in the framework of microscopic theory in a good 
agreement with experimental datatl. However, it does not 
elucidate the nature of the processes resulting in the re- 
construction of the structure at the melting point. There 
are direct evidences of this reconstruction such as nonlin- 
ear in temperature T increase of the heat capacity and a 
significant softening of shear moduli at T — > Tm (for the 
most of metals the latter is as large as 50% decreasej-pf 
one of the moduli but, e.g., for indium it is even 90%Q). 
To investigate the phenomena taking place in the crys- 
tal lattices near T^ a broad variety of techniques are 
used - froin-|dassical or quantum Monte-Carlo computer 
simulationslZlapto experiments with artificial "crystals" of 
charged dropaa. Despite this we have no yet a complete 
and reliable picture of the melting of crystals. Here we 
try to suppose a qualitative scenario of the melting for a 
specific case of metals. 

We are interested in the dynamic picture of the pro- 
cesses in a crystal near the melting point. According to 
the well-known Lindeman criterion of the melting one 
has u/d ~ 0.1 at T = Tm where u = ■^/(u^), (v?) being 
average square of thermal atomic-displacements, d is the 
minimal interatomic distance (seed) . It is important that 
under such conditions, as a rule, phonons are well-defined 
collective excitations in all the Brillouin zone and anhar- 
monic effects are relatively small up to Tm- Microscouk: 
calculations for metals with BCC and FCC structuresE2l 
show that the value of relative phonon damping Fq^/wq^ 
(where q,^ are the wavevector and phonon branch num- 
ber, correspondingly) at T = Tm does not exceed 0.1-0.25 
for the most of them. It is naturally to ask why so small 



anharmonic effects may lead, nevertheless, to the loss of 
stability of crystal lattices at the melting. To our opin- 
ion, the answer is connected with geometric properties of 
three-dimensional Euclidean space. 

Microscopic calculations of structural apd thermody- 
namic jxharacteristics of liquid rare gasesllil and liquid 
metalgj show that the accuracy of their description near 
Tm is dependent mainly on the simple geometric param- 
eter of atomic packing ry = T:d^n/& where n is the atomic 
density, d is the diameter of non-overlapped spheres con- 
nected with each atom. It confirms qualitatively the 
adequacy of geometric approach to the_xlescription of 
structure of liquids proposed by Bernallj, namely, the 
model of random packing of hard spheres. His conclusion 
about the existence of close packed "Bernal pseudom: 
clei" formed by the connection of several tetrahedraliScI' 
is the most important for the scenario of the melting 
proposed here (they are called pseudonuclei because it 
is impossible to build regular structure from these for- 
mations). Locally they are characterized by the packing 
parameter 77 ~ 0.78 (at the average packing parameter 
77 ~ 0.64) which is higher than the maximum possible 
packing in crystals, r] ~ 0.74 (seeu). To prevent mis- 
understanding note that we consider purely geometric 
definition of the packing parameter corresponding to the 
hard-sphere model with the interatomic potential 

Such definition does not take into account the effects of 
"softness" of interatomic interactions in metals, i.e. the 
role of the attractive part of the interactions, as well as 
nonpairwise forces (the dependence of the potential on 
mean electron density) which are important for metals. 
Because of this the effective packing parameter in liquid 
metaJs at T ~ Tm can be estimated really as rjmet — 
0.45lI. Despite this molecular dynamics simulations give 
us evidences about qualitative adequacy c^ the Bernal 
model of the structure of liquids as a wholeB. 

For metals as well as for liquid rare gases the clos- 
est packing is energetically favorable for small groups of 
atoms what is confirmed byrthe data about the struc- 
ture of small metallic clusteralj. It is well known that in 
three-dimensional Euclidean space the closest packing is 
the tetrahedral one which, at the same time, cannot be 
reached in all the space. This is a drastic difference of 
three-dimensional space from the two-dimensional one. 
The latter can be filled completely by regular triangles, 
six of them being meet in each lattice site. In three- 
dimensional case only five regular tetrahedra can have 



the common edge and the void appears in this case with 
the angle deficit S = 2n-5 cos"^ (1/3) , d/2iT ~ 0.02. This 
simple consideration is the base of an elegant approach 
to the description of structure of liquids and glassesQ. 
According to this approach the structure is defined not 
with the respect to usual Euclidean space but to Rieman- 
nian one. The latter can be, for a specific value of the 
curvature radius, filled completely by regular "tetrahe- 
dra" . The transition to the real physical space is carried 
out by introducing of corresponding "decurving" defects, 
namely, structural disclinations. Their appearance near 
Tm has been recently demonstrated in a simple, model 
of a crystal by molecular dynamics simulationsEil. Ear- 
lier a similar approach to the problem of melting bas- 
ing on the consideration of statistics of linear defects has 
been developed phenomenologically by Patashinskii with 
collaboratorst£l 

Taking into account the elastic energy of the disclina- 
tions as well as the trends to the closest packing of atoms 
in "chemical" energy of interatomic interactions it may 
be showntEl that the formation of an inhomogeneous state 
can be in principle described by these considerations. 
This state is characterized by the existence of super- 
dense packed regions and Bernal "voids" around them. 
The central regions have tetrahedral-like packing struc- 
ture and are similar in this sense to Bernal pseudonuclei. 
Our main hypothesis is that such regions can be created 
by thermal fiuctuations near the melting point. 

To form the close packed region as a result of thermal 
motion of an atomic group it is necessary to "overcame" 
the angle deficit 6. It has to be done as a result of not 
just oscillations but of displacements of the centers of the 
oscillations, i.e. of long-lived component of atomic mo- 
tion. The latter can be estimated from the following sim- 
ple considerations. The irreversibility of atomic displace- 
ments results from the phonon damping. The weight of 
this irreversible component is of order of T/uJ where F, U 
are the average phonon damping and frequency, corre- 
spondingly. The close packed region will be long-lived 
enough after its formation by thermal fiuctuations since 
it corresponds to a local minimum of free energy. In- 
deed, it is known that for a system of hard spheres with 
relatively weak attraction |-(|which is a good approxima- 
tion for the most of metalsErEIl) the density of free energy 
of sufficiently small groups-of. atoms is minimal for the 
closest tetrahedral packin^ll3. As a result, the proba- 
bility of formation of the close packed pseudonucleus can 
be not too small for atomic displacements satisfying the 
condition 

{u/d) (r/ZU) ~ 6/27r ~ 0.02. 

We believe that this estimation is a realistic melting cri- 
terion. Empirically, it coincides approximately with the 
Lindeman criterion but, in contrast with the latter, it is 
directly connected with contemporary views about the 
structure of liquids and a scale of anharmonic effects, at 
least, for metals. 



Provided that the average atomic density is fixed the 
formation of close packed region results inevitably in the 
appearance of Bernal voids around it. This changes es- 
sentially traditional views formulated by FrenkeB about 
the leading role of the monovacancies in the formation of 
free volume near the melting point. According to the sce- 
nario proposed the latter appears first with the temper- 
ature increase as the voids around close packed regions. 
Molecular dynamics simulation shows that for systems 
with "soft" enough interatomic interactions (e.g., liquid 
metals) there are two kinds otiBernal voids: distorted 
tetrahedral and octahedral oneso. Apparently these voids 
are the most energetically profitable "carriers" of the free 
volume in the solid phase also. Monovacancies should 
have higher free energy and appears only as a result of 
diffusion decay of the voids in the close vicinity of Tm ■ 

To check this assumption we have made the fitting 
of defect contributions to the heat capacity of sodium 
near Tm which can be obtained from the experimen- 
tal data by substracidon electronic, phonon and anhar- 
monic contributionalj. Trying to fit this contribution 
C^ (T) by Arrhcnius law in a broad temperature range 
(0.6Tm < T < T„i) we find that the decline of the curve 
has a jump at T 2± O.STm- The fitting for the narrower 
temperature range O.ST™ < T < Tm gives the follow- 
ing results for the entropy Si^ and energy Eiy of the 
formation of monovacancies as well as their concentra- 
tion at the melting point ni^: Siv — 4.3, -Biu ~ 0.3 eV, 
niy ~ 0.3%. For another temperature range O.GTm < Tj-< 
0.8r,„ one has Siy ~ 0.5, Si„ ~ 0.17 eV, ni„ ~ 0.74%EI. 
The first set of the results is in much better agreement 
with the data obtained by differential dilatometry data 
Sly ~ 5.2, Ely ~ 0.35 eV, ni„ ~ 0.3VM as wefi as the 
results of microscopic calculations of |tiie energy of for- 
mation of monovacancy, Eiy ~ 0.3 eVo. 

This confirms that, in contrast with traditional views 
and in agreement with the scenario of the melting pro- 
posed here, the temperature increase leads first to the 
formation of clusters of vacancies (the voids) and the 
monovacancies appears only in the close vicinity of Tm- 
Note also that the formation of monovacancies in a broad 
temperature range below Tm corresponding to Frenkel 
views about the nature of the melting was never oh^ 
served in computer simulations by molecular dynamicsti 
and Monte-CarloQ methods. At the same time, according 
to these simulations, pentagonal faces of Voronoi polyhe- 
dra are observed frequently for instant atomic configura- 
tions notjOHly in liquids but also in crystals near the melt- 
ing pointcHZl. According to Bernal this is the main feature 
of the structure of liquids. It means that some struc- 
tural elements which are characteristic for liquids arises 
by thermal fluctuations already in the crystal phase. 

As a conclusion, consider a possible way to describe 
such fiuctuations quantitatively. The formation of the 
closest tetrahedral packing near Tm has toJead to the 
appearance of icosahedral short-range ordeiu describing 
by the parameter 



Q6m(r)=F6m[0(r),^(r)] 

where Ygrn are spherical functions, 9 (r) ,if (r) are po- 
lar angles of the direction to the nearest neighbors of 
the atoip-,at point r and the corresponding correlation 
functionEJ 

GW = SE(Q6'»WQ6™(0))- 

m 

Its calculation in the crystal phase by molecular dynam- 
ics or Monte-Carlo methods could give an important in- 
formation about the local structure of solids near the 
melting. 

This work is supported by Russian Basic Research 
Foundation, grant 98-02-16219. 



(1989) [Engl. transL: Sov. Phys.- Solid State 31, 427 

(1984)] 
^'^ V. G. Vaks, S. P. Kravchuk, and A. V. Trefilov, J. Phys. 

FIO, 2325 (1980); Uspekhi Fiz. Nauk 155, 159 (1988) 

[Engl. transL: Sov. Phys. -Uspekhi 31, 474 (1988)] 
^^ W. Adihart, G. Fritsch, and E. Luescher, J. Phys. Chem. 

Sol. 36, 1405 (1975) 
^^ A. M. Bratkovsky and N. E. Zein, Fizika tverd. tela 26, 

2561 (1984) [Engl. transL: Sov. Phys.- Sohd State 26, 1553 

(1984)] 
^^ P. J. Steinhardt, D. R. Nelson, and M. Ronchetti, Phys. 

Rev. B28, 784 (1983) 



^ J. Frenkel, Kinetic Theory of Liquids (Clarendon Press, 

Oxford, 1946) 
■^ V. P. Skripov and V. P. Koverda, Spontaneous Crystalliza- 
tion of Supercooled Liquids (Nauka, Moscow, 1984) 
^ S. M. Stishov, Uspekhi Fiz. Nauk 154, 93 (1988) [Engl. 

transL: Sov. Phys. -Uspekhi 31, 52 (1988)] 
* M. Kleman and J. F. Sadoc, J. Physique- Lett. 40, L569 
(1979); D. R. Nelson, Phys. Rev. B28, 5515 (1983); M. 
Kleman, Adv. Phys. 38, 605 (1989) 
^ A. M. Bratkovsky, V. G. Vaks, and A. V. Trefilov, ZhETF 
86, 2141 (1984) [Engl. transL: Sov. Phys.-JETP 59, 1245 
(1984)] 
^ V. E. Zinoviev, Properties of Metals at High Temperatures 

(Metallurgia, Moscow, 1989) 
'■ L. Gomez, A. Dobry, and H. T. Diep, Phys. Rev. B55, 

6265 (1997) 
® M. Boninsegni, C. Pierleoni, and D. M. Ceperley, Phys. 

Rev. Lett. 72, 1854 (1994) 
''H. M. Thomas and G. E. Morfil, Nature 379, 806 (1993); 
A. E. Larsen and D. G. Grier, Phys. Rev. Lett. 76, 3862 
(1996) 
^° V. G. Vaks, S. P. Kravchuk, and A. V. Trefilov, J. Phys. 
FIO, 2105 (1980); M. I. Katsnelson, A. V. Trefilov, and K. 
Yu. Khromov, Pis'ma ZhETF 69, 649 (1999) [EngL transL: 
JETP Lett. 69, 688 (1999)] 
" J. A. Barker and D. Henderson, Rev. Mod. Phys. 48, 588 

(1976) 
^2 J. D. Bernal, Nature 183, 141 (1959) 

" L. D. Marks, Philos. Mag. A49, 81 (1984); M. W. Sung, R. 

Kawai, and J. H. Weare, Phys. Rev. Lett. 73, 3552 (1994) 

^^ P. Jund, D. Caprion, J. F. Sadoc, and R. JuUien, J. Phys. 

: Cond. Matter 9, 4051 (1997) 
^^ A. C. Mitus and A. Z. Patashinskii, Physica A150, 371, 
383 (1988); A. Z. Patashinskii and L. D. Son, ZhETF 103, 
1087 (1993) [EngL transL: JETP 76, 534 (1993)] 
^^ Yu. N. Gornostyrev, M. I. Katsnelson, and A. V. Trefilov, 

J. Phys. : Cond. Matter 9, 7837 (1997) 
^■'V. G. Vaks, V. G. Kapinos, Yu. N. Osetskii, G. D. 
Samoliuk, and A. V. Trefilov, Fizika tverd. tela 31, 139 



